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- We assume that the properties of nonequihbrium stationary states of systems of particles can 

r> I be expressed in terms of weighted orbital measures, i.e. through periodic orbit expansions. 

' This ahows us to derive the Onsager relations for systems of particles subject to a Gaussian 

thermostat, under the assumption that the entropy production rate is equal to the phase space 
contraction rate. Moreover, this also allows us to prove that the relevant transport coefficients 
are not negative. In the appendix we give an argument for the proper way of treating grazing 
collisions, a source of possible singularities in the dynamics. 



1 Introduction 

Very recently a number of proofs have been presented of the Onsager reciprocal relations (OR) 
in non-equilibrium stationary states based on dynamical systems theory [|l|, |2|, ^. There, the 
earlier assumption of Refs.|Q, ||], i.e. that the dynamics of systems of particles could be re- 
garded as "Anosov" for practical purposes,tJ was used to characterize the relevant stationary 
distributions. Consequently, such distributions were conjectured to be of the kind first investi- 
gated by Ya. Sinai, D. Ruelle and R. Bowen (SRB distributions) [^. This assumption makes 

^For completeness, the main assumption of Ref.|^, p. 935 states tliat: 
A reversible, many particle system in a stationary state, can be regarded as a transitive Anosov system, for the 
purpose of computing its macroscopic properties. 
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possible a generalization to nonequilibrium stationary states of the classical equilibrium ensem- 
bles for the probability to find a system in a certain phase in phase space. In particular, this 
can be done by attributing to the absolute values of the Jacobian of the dynamics restricted to 
the unstable manifolds, a meaning connected to that of the Boltzmann factor of the canonical 
ensemble (see subsection 2.2 below, and see Ref.|0] for numerical experiments meant to justify 
the connection). Use of such techniques lead to a fluctuation theorem —first investigated in 
Ref.|^, ^, m and further confirmed by a variety of computer experiments [P, IC]— and to the 
above mentioned proofs of the OR. 

There are two representations of the SRB distribution which have been used in practice 
so far: the Markov partition method of Ref.|||], and the periodic orbit (or cycle) expansions 
(POE) (see, e.g. Refs.fl^, |l2|, ||]) used, for instance, in Refs.^, |l6|, [l^, [11. Here, we 
present an approach to the proof of the OR based on the assumption that the stationary 
distribution can be expressed in terms of POE, and on analyticity conditions on the properties 
of the relevant orbits. The point is not to improve on the quite general proof of Ref.^], but 
to show that a different approach with different assumptions can lead to similar results, thus 
shedding new light on the precise connection between the OR in statistical mechanics and the 
dynamics of mechanical systems. Periodic orbit expansions are also used, here, for a proof 
that the transport coefficients of our systems are nonnegative, following in part ideas developed 
previously in Ref.|jl^ (see also the paper by D. Ruelle, Ref . [pO|] ) based on the equality of phase 
space contraction rate and entropy production rate. 

The possibility of describing the stationary states of systems of particles in terms of POE 
is implied by the assumptions of Refs.Q, |5|,0 but the actual range of applicability could be 
different (see next section). The importance of this approach lies in the fact that, although 
the practical implementation of the POE for a system of many particles is exceedingly difficult 
(practically impossible with present day technology), its mathematical expression can be used 
for formal proofs of properties of real systems. Moreover, the validity of given results obtained 
in terms of POE can be tested directly on simple models of statistical mechanical interest, 
by means of modern computers. In fact, POE have been succesfully applied in numerical 
simulations of one particle systems, such as various (periodic) Lorentz gas models, see e.g. 
Refs.[0, 15, [l^. In this respect, it is also important to note that the efficiency of this 

kind of numerical studies can be greatly enhanced, as recently shown in Ref.[21|, thus making 
the numerical implemetation of POE conceivable for more complex systems. 



2 Properties of transport coefficients from POE 
2.1 Preliminaries 

In this paper, we concentrate on a d-dimensional system of N interacting particles, subject to 
external fields of force (whose strength is gauged by parameters), and subject to a Gaussian 
thermostat |2^. Moreover, in order to make use of the techniques of dynamical system theory, 

In other words, if the (possible) deviations of the dynamics from the Anosov conditions bear no consequence 
for the calculations of the quantities of physical interest, they can be considered as irrelevant. In that case, the 
system will be called "Anosov-like" . 

^If the system under study is really Anosov, weighted averages of smooth functions evaluated over periodic 
orbits do converge to the averages computed with the SRB distribution (see e.g. Refs.||ll|, |l^, Q), despite 
the fact that periodic orbits cover but a set of zero (Lebesgue and natural) measure. If the system is only 
"Anosov-like" , POE should work as well, at least for the calculations of the physically interesting quantites. 
However, in this paper we will assume the validity of POE only for the calculation of the average currents, in 
order to limit the assumptions as much as possible. 
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we impose periodic boundary conditions on our models. In particular, the equations of motion 
we consider take the form 

qi = — , Pi = *»(q)+Fi(G)-a(p,q)p, , i = l,...,iV, (1) 

where a(p, q) characterizes the Gaussian thermostat and is determined by the constraint one 
wants to implement (e.g. fixing the value of a special dynamical quantity). Also, qj, pj and 
rrii are respectively the coordinates, the momentum and the mass of the ith particle, with 
(P)q) = ((Pi)i^i) (qi)i^i); is the force on particle i due to the interactions with the other 
particles; and Fj is the action of the external fields on the ith particle, which depends on 
the values of the parameters G = (Gi, ...,0^)- The purpose of the Gaussian thermostat is to 
make such an A^-particle system reach a stationary state when the fields are "switched on", 
starting from a given initial distribution in phase space. We assume that when G vanishes, 
the external forces Fj and the thermostat coupling a also vanish, making the stationary state 
an equilibrium state. 

The particular constraint we impose on our models is that the total energy remains fixed 
in time (i.e. we use an iso-energetic thermostat). This yields 

«(r;G) = 4.E^^^^^; i^ = E— (2) 

1=1 ' j=i ' 

where F = (p, q) stands for a generic point in the phase space, and K is the total kinetic 
energy of the system. S Once the stationary state has been achieved, the periodic boundary 
conditional imply that there is an elementary cell (EC), O say, whose replicas cover the whole 
(2(iiV)-dimensional phase space, and there is a time evolution on 17, denoted by St, which 
represents the dynamics of the N particles in 0. One may also consider the dynamics with 
respect to a given Poincare section, V say, and the timing events for the definition of V may be 
chosen to be the collisions between particles. When the interacting potentials are soft, this can 
be done by calling a collision the event in which two or more particles come within a certain 
distance of each other |^] . 

By identifying the phase space contraction rate for the dynamics —the divergence of the 
right hand side of Eqs.([l|)— with the (microscopic) entropy production rate, we define the v 
(microscopic) currents as 

J(r; G) = VGa(r; G) ^ (Ji, J,)(r; G) . (3) 

This definition of current was first adopted in Ref.[^, and also used in Refs.|^, |l9| ]. 

In order to develop our analysis in terms of POE, we need that Unstable Periodic Orbits 
(UPOs) be densely embedded in the support of the natural measure of the system and, if they 
are,0 we must characterize their properties. In the first place, UPOs are orbits with at least one 
positive Lyapunov exponent, and it is convenient to collect them in sets of orbits with a fixed 
number of collisions. For instance, denote by Pn{G) the set of orbits within which n collisions 
occur. Then, we let uj{G) indicate a generic element of Pn{G), and Tj^(g) its period. Given 

^In case of hard core particles, the iso-energetic constraint coincides with the iso-kinetic one, i.e. with the 
constraint of constant kinetic energy. 

^Periodic boundaries may look as a rather restricitve condition. However, one should keep in mind that Vl 
can be arbitrarily large and complex, as long as it may be used to periodically tile the infinite phase space. 

^Typical dynamical systems whose UPOs are dense in the support of the natural measure are those verifying 
the ajdom-A conditions. Indeed, for such systems the density condition is part of their very definition. 
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(a) (b) 

Figure 1: The equivalence between the dynamics in the EC and those in the full phase space, 
a) Trajectories of two black round particles in the infinite phase space where they move never 
returning to the same point. Only two adjacent cells, out of the infinitely many which tile the 
phase space, are drawn, b) Representation of the situation (a) in the EC, where the orbits 
appear to be periodic. To average any quantity over the trajectories in (a) is the same as 
to average that quantity over the trajectories in (b), as long as the actual diplacements are 
accounted for. 



any orbit cj(G), with Lyapunov exponents Xuj{g),i ^ Xu)(G),2 ^ ••• ^ '^c<;(G),2d7V) we consider the 
following quantity: 

K{G),u = exp {-Tu^iG) E ' ^^{G),j) > (4) 

where Xuj{G),j is the j-th Lyapunov exponent of orbit and the prime on the symbol J2 

indicates that the sum runs over the positive exponents only. The subscript u on the left 
hand side of Eq.(^) says that A.^(^g),u is computed restricting the dynamics to the unstable 
manifold. Notice that A^(^g),u is a measure of the instability in the phase space of uj{G), and 
it constitutes the basic ingredient in any POE, i.e. it is part of the (unnormalized) weight 
attributed to uj{G). In particular, the less unstable the UPO, the larger its weight, as explained 
in Section 2.2. 

For the study of transport in our particle systems, it is interesting to consider the total 
displacement of the single particles during a whole period of a UPO. In this respect, given 
io{G), it is important to recall that the displacement of the jth particle at the end of one period 
of uj{G), Aqj ,^((3), does not necessarily vanish. When it does not, also the average current 
associated with uj{G) may not vanish. Here the average current is the integral of Vqck over 
the UPO (that is, the cumulative current I^^ defined below in Eq.(|5|)), divided by the period 
Ti^(G)- -^or t'^is reason, it makes sense to talk of currents associated with a UPO although, at 
first sight, it may look strange. Indeed, the periodic boundary conditions impose translation 
invariance symmetry along the directions of the lattice vectors in the stationary state. Hence, 
to follow a trajectory which never returns to the cell where it started from, is equivalent to 
following its translated image in the EC. See, e.g., Fig.l: the trajectories of the round black 
particles in the full phase space do not close. However, they can be followed in the EC, where 
they appear to close producing a periodic orbit. Clearly, to average quantities over the open 
paths —the real paths in phase space— is the same as averaging them over the images of such 
paths in the EC, as long as one keeps track of the actual distance covered during the motion. 
In this sense, when we speak of periodic orbits, we mean all those that are such in the EC. In 
particular, if the particles in Fig.l carry the same charge c, the side of the EC is i, and the 
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period of the orbit is r, the relevant current is 2c£/r. We remark again that it is the presence 
of the periodic boundary conditions which, in the stationary state, ahows us to map the infinite 
system into a finite one: the EC. Away from the stationary state, or in the case that we do 
not have periodic boundary conditions (hke in the random Lorentz gas), the techniques of 
dynamical system theory fail to apply in this simple manner. 



2.2 Conditions 

The present paper rests on Assumptions 1. and 2. below which, similarly to those of Refs.|Q, ^, 
are empirically motivated. In particular, our assumptions were inspired by the numerical 
results presented in Refs.|0, 15, 17, |l^] and, in a broader sense, by those of Ref. [p. 

Let us introduce the quantity laj(G)) the cumulative current associated with uj{G) G 
Pn{G): 



I^(G) 



/ J(5tr,(G)) dt , (5) 

Juj(G) 



where the integral is carried over one period of uj{G), and StT^(^Q^ is the point representing 
the state of the system at time t, if it was ^uj{g) ^ ^{G) at t = 0. Then, our first assumption 
can be formulated as follows. 

Assumption 1 (Existence and representation of stationary states). The dynamical 
system {Vt, St) describing a reversible many-particle system obeys the "Extended zero-th law" 
of Refs.J^ at least for what concerns the currents. Then, letting fic the corresponding 
(SRB) stationary distribution, we can write 

(J)„(G) = hm i f^JiStT) dt= [ j(r) df,G{r) , (6) 

T^oo I Jo JQ 

for (Lebesgue) almost all T gVL. Also, the phase space average of J with distribution fic can 
be given in terms of the following POE: 

(J),(G) = hm ^^^^.iWr. ^ ^.^ E.gP4G) ^^!ul^ 

Equation (|^ expresses the current as a limit of weighted averages of orbital average currents, 
where the weights have the suggestive form 

r^Aj^ = (time spent in oj) x (inverse of instability of oj) , (8) 

apart from a normalization factor. This weight associates the points of a given region of phase 
space with a probability and is larger for longer UPOs, while it is smaller for more unstable 
UPOs. It is in this sense that Tf^A"^ plays a role similar to that of the Boltzmann factor in 
the canonical ensemble. 

Let us denote by (J)^"^ the "n-th order UPO approximation" of the current, i.e.: 



(J)(")(G) = ^"^''"^^)'^"-"_^" , (9) 



where all properties of a given UPO, uj say, depend on G, because oj is taken from the set 

Pn{G). 
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Assumption 2 (Properties near equilibrium). There is a positive constant p such that 
for [G| (the modulus of G) smaller than p we have: 

a. the period t^(q), the stability weight -^^c) u' ^^'^ thermostat coupling a{StT^; G), are 
"left" and "right" analytiM with respect to all the components ofG, for alluj{G) € Pn{G), 
all n € IN , and all t > 0. 

b. (J)*-"-' (G)}, converges uniformly to a limit as n ^ oo, for i = 1, 

c. The support of p,G is dense in 0,. 

In particular, this assumption implies the existence of constants a+, aj, b+, b^^, c+ and in 
IR" such that 

T-c(G) = Toj{0) + • G + o{G'^) , (10) 

KtG),u = Kku + ^^-G + o{G'), (11) 

»{Str^iG)) = c^miStT^m) ■ G + o(G2) , (12) 

where the i-th component of r^^, r^^i say, for r = a, b, c, equals r^j when Gi > 0, and r^^i = rj • 
for Gi < 0. The necessity for left and right analyticity emerges in the case of special geometries, 
i.e. as a consequence of the actual shape of the particles and of the medium where they move. 
There are cases, indeed, such that reversing the sign of the components of G may require a 
change in the values of the constants a, b and c above, so that only the left and the right 
analyticity around Gi = 0, for all i, hold. For our purpose, this is sufficient, but in the case 
we have full analyticity at G = 0, some approximations can be improved, as explained below, 
after Theorem 1. 

Assumptions 1. and 2. are needed here to justify the calculations of subsections 2.3 and 
2.4 but, similarly to the assumptions of Refs.Q, |5|, they cannot be validated at present on the 
sole grounds of the given dynamics. Only a posteriori we can check whether using them we 
are led into some kind of inconsistency or not. However, to accept them as possibly valid (a 
priori) we may rely on the evidence accumulated so far in the literature, e.g. on several studies 
of the nonequilibrium Lorentz gas like Refs.|16, 23| and others. In particular. Assumption 



2. a is also in agreement with the discussion of pp. 949-950 of Ref. [g], where it is argued that 
the transition from equilibrium to nonequilibrium stationary states, with small forcing, should 
be seen as producing "... an insignificant deformation of the unstable manifold Wq i.e. of 
the unstable manifold of a given periodic point. Similarly, one may argue that the properties 
of periodic orbits will only change little for small changes in the fields. As a matter of fact, the 
periodic (trianglar lattice) Lorentz gas has been used to test many orbits to high numerical 
precision in order to see whether their contributions to numerator and denominator of Eq. 



are smooth in the field. The result is that they appear to be smooth |18|, with the only 
exceptions of orbits with grazing collisions, i.e. orbits which suddenly appear or disappear 
from the phase space, when the field is varied ever so slightly. However, it does not seem 
that such orbits should concern us. In the first place, they do not exist if the interactions are 



^Given a function f : M ^ ]R, and a constant p > 0, we say that / is left analytic in {—p,0] if there are 
constants /o and r~ in K such that f{x) = fo + r~x + o(x^), for all s G (— p, 0]. Similarly, we say that / is right 
analytic in [0, p) if there are constants /o and in R such that f{x) = fo+r^x + o{x^), for all x £ [0, p). Thus, 
a function which is analytic in a neighborhood of zero is left and right analytic, and has r~ = r^. However, / 
could be left and right analytic without its derivatives being defined at a; = 0. Ifere o{x^) is to be understood 
as a term of second order in x, negligible with respect to the other terms for all x in the relevant intervals. 
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soft, and hard core interactions could be seen as a limiting case of soft core ones, Ref.|| 

Also, the numerical tests of Refs.0, 17, 18 1 have not evidenced any difficulty linked to this 



problem, because the orbits with grazing collisions were a very small fraction of the whole, 
with negligible weight each. In the Appendix, we give an argument to justify why these orbits 
should always be assigned a vanishing weight in the expansion Eq.(0). 
A few further remarks are in order. 

Assumption 2. plays for us the role of the differentiability of SRB measures in Ref.|p, but 
it may appear as rather strong. On the other hand, this assumption may be relaxed in various 
ways, and it is supported to a good extent both by numerical studies of simple systems, and 
by intuition, as noted above. In practice, we observed that for small changes in the field the 
orbits change very little in shape which, in turn, determines their properties such as their 
period, stability, a etc. Furthermore, it is important to note that we do not require that our 
conditions hold for very general functions of phase. On the contrary, we restrict the validity of 
our conditions to very special quantities —the phase space contraction rate and its derivatives— 
which are easily seen to have particularly good properties in the dynamics of many interesting 
maps. 

A comparison between Assumption 1. and the assumptions of Ref.|Q shows that the latter 
directly refer to the dynamics of the system, i.e. to the properties of the equations of motion 
and of the space in which the motion takes place. Our assumption, instead, refers to the 
stationary measure, which is not the dynamics but only a result of the dynamics. Therefore, 
Assumption 1. is valid if the system is Anosov or axiom-A (hence, also if the assumptions of 
Ref.||5| hold) but it could still remain valid even for dynamics of a different kind. On the other 
hand, our assumption is not as fundamental because, for instance, it does not provide us with 
a direct way of estimating the errors connected with the n-th order UPO approximations of 
the POE. 

2.3 Derivation of the Onsager reciprocity relations 

We can now prove the validity of Onsager's relations for the systems verifying our assumptions 
(which we take to be general A^-particle systems subject to an iso-energetic Gaussian thermo- 
stat). The proof proceeds as follows. First, we consider the n-th order UPO approximation of 
Eq.(0), and we expand to first order in G the resulting expressions. Then, we group together 
the contributions of periodic orbits with opposite currents, in order to use the relation called 
Lyapunov sum rule in Ref.|l|]. This is a relation between the parameter a and the Lyapunov 
exponents associated with a given periodic orbit (see Eq.(^) below). This allows us to find a 
symmetry between different derivatives of the n-th order UPO approximate current. Finally, 
exploiting Assumption 2. we interchange the derivative with taking the infinite period limit 
operation, to obtain the desired result. The major ingredient in these calculations is the time 
reversal symmetry of the equations of motion, together with the density of the attractor in 
the phase space. Indeed, it is because of these two features that we can write Eqs.(|l|,|l5|,0) 
below. 

Using the assumed analyticity in the fields of the period and of the stability weight 
(Eqs.(|lO|,|lT])), we can expand the denominator of Eq.(|9|) to first order in G and obtain 
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(J)W(G) = 



E, 



E 



G + o(G2 



(13) 



Then, the time reversibihty of Eqs.(|T|), combined with the density of the relevant attractor in 
phase space, guarantee that for every orbit co £ EpiG) there is another orbit —lo £ Pn{G) with 

1,0 i.e.: 



equal period and opposite cumulative current 



-u){G) — Tui(G) ; 



Hence, by grouping together such pairs of orbits, we get: 



^eP„(G) 



-1 



E 

a..GP+(G) 



A, 



-1 



A 



-1 



E 

L.6P+(G) 



IljAj^ ,^ (1 



(1-A^) , 



(14) 



(15) 



where uj G P^{G) if it has n collisions and a{StT^)dt > 0,1 and 



A, 



^ = exp 



2dN 
E 



(16) 



where the sum this time involves all the Lyapunov exponents of uj. The Lyapunov sum rule 
for a periodic orbit to -Eq.(8) in Ref.|l£]- in our context can be written as 

2dN 



-[dN -I) 
-{dN 



u;(G) 



1) / c^^o){Str^{o))dt-G + o{G' 

Ju){G) 



(17) 



Note that for G 
A, 



the flow does not contract or expand elements of phase space, i.e. 
POs, and each addend in the sums in Eq.(|l5|) vanish. Then, expanding 
the terms in Eq.(|l5|) in powers of G, and substituting into Eq.(|l3|) we obtain 

(J)(")(G) = A(")G + o(G2) , (18) 

where A^") is the symmetric matrix whose ij-entry is 

Et^GP„(G) ^t^('o),«Qt^{0),jQa;{0)j 



A. 



in) 



{dN - 1)- 



E 



ujePniG) ■^t<^(o)A^(o),« 



(19) 



^Observe that time reversal symmetry in the equations of motion is not enough for our purposes (see e.g. 
Ref. [^5|). The equations of motion of our models are time reversible for all G, however the corresponding 
attractors and repellers may cover disjoint regions of phase space if G is sufficiently large. In that case, taking 
the time reverse image of a point in the attractor produces a point in the repeller, and it is not true then that 
—u! G Pn{G) if oj G P„{G). On the contrary, the time reverse of a UPO in Pn{G) is still in P„{G), if the 
attractor is dense (as in Anosov systems). Systems with non dense attractors can be Studied (seee.g.Refii), 
but they do not concern us here, as we are interested in the small fields regime. 

**Given n G iV, it may happen that there are no UPOs in Pn{G) with ^^a{StVi^)dt > 0, in which case 
P^{G) is empty and the sums in Eqs.([l^ vanish. If this happens for all n greater or equal to a given no, then 
{J)n = 0, and our calculations leading to Eq.(E3|) trivially hold. 
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and (5cj(o),i is defined by 

Qw(o),i = / Cj^(o),j('S'tr^(o)) dt . (20) 

Tlius, if we differentiate the A;-th component of the approximate current, (Jfc)^")(G), with 
respect to Gi say, we obtain 

aGXJ.)(")(G) =4?+oi(G) =4")+02(G) =9g,(J^)(")(G)+03(G) , (21) 
where the Oi{G), i = 1,2, 3, are distinct quantities of order G. Let us denote by 

Lm = dG,{Jk)n\G=o = dc, ^^^Jh)^\^^ ' (22) 

the /cZ-entry in the transport coefficients tensor, and recall that the sequence {{Jk)^'^\G)}'^^i 
converges to a limit (Assumption 1) uniformly (Assumption 2), as n — > oo. Then, we have 
(Assumptions 1. and 2.) 

Lki = da, lim ( Jfc) W ^ = hm da, ( Jfc)(") =Lik\ k,l = 1, u . (23) 



G=0 



These are the Onsager relations for our systems, and we can state that: 

Theorem 1 For thermostatted N particle systems verifying Assumptions 1. and 2. the fol- 
lowing holds 

Lki = Lik for all k,l = 1, d . (24) 

As mentioned above, we note that the geometry at hand may improve our calculations, 
thanks to ensuing symmetries of the corresponding systems. By geometry we mean all the 
specifications which have to do with the shape of the particles and of the medium in which 
they move. In most cases, these are such that reversing the direction of G effectively amounts 
to just a rotation of the coordinate axes. In other words, the UPOs embedded in the attractor 
corresponding to — G have the same properties as those in the attractor corresponding to 
G, except for having an opposite cumulative current: i.e. for every n E IN, and for every 
uj{G) e Pn{G) there is uj{-G) G P„(-G) such that 

M-G) = mg) ; K(-G),u = K(G),u ; M-g) = -Mg) ■ (25) 

Observe that this has nothing to do with time reversal invariance, although it may appear to 
be similar. Equations ( [25| ) are a consequence of the geometry of the system only, and express 
what can be referred to as the [G — > — G]-field reversal symmetry, which may not obtain for 
anisotropic systems. Therefore, having distinguished the time reversal invariance from the 
(possibly not present) field reversal symmetry, we see that the second is not needed for the 
validity of the Onsager reciprocity relations. If, on the other hand, Eqs.(p5[) hold. Assumption 
2. can be taken in the sense of full analyticity,0 and the terms I^(g')(A^^^ ^ — A~^j.gj ^) of 

Eq.([l5|) are odd in G, while the terms Tuj(g)^ZIg)u -^^-(H) ^^^^ even in G. Then, Eq.(|l8|) 
can be re-written as 

(J)(") (G) = A(")G + o(G3) (26) 

making our calculations correct to second order, rather than just first order in the field. 
However, such field reversal symmetry is not needed for the derivation of the OR. 



'For instance, in Ref. full analyticity of a around G = is assumed. 
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2.4 Nonnegativity of transport coefficients 



Here, we generalize an argument recently proposed in Ref.[19|. Let us consider the transport 
coefficient associated with the k-th current: 



cTk = dGk{Jk)n\G=o = Lkk , k = l,...,u. (27) 

Then, in view of Eq.(^), and grouping terms as in Eq.(15), allows us also to deduce the fact 
that CTfc cannot be negative. To see that, rewrite Eq.(0) as 

(Jfc)„(G)=lim ^ — ^ = Y,AkiGi + o{G') ] k = l,...,u (28) 

which implies ak = A^k- This quantity, in turn, is nonnegative as it results from taking the 
limit A^k = liiHn— »oo 

A^^^ in Eq.dll) (the limit exists because of Assumptions 1. and 2.). 
Therefore, we can state this as 

Theorem 2 For thermostatted N particle systems satisfying Assumptions 1. and 2, > 0, 
for all A; = 1, v. 

For this result Assumptions 1. and 2. are enough. On the contrary, the problem of the 
strict positivity of does depend on further conditions, similarly to the situation investigated 
in Ref.pl. 



We conclude this section by observing that all our results are not purely dynamical in 
nature: they rest on the choice of the initial ensemble, which must evolve into the stationary 
SRB measure by the given dynamics. Indeed, initial conditions such that the consequent 
evolution violates our conclusions are possible. However, Assumption 1. says that such initial 
conditions are very special, because they only constitute a set of zero Lebesgue measure. We 
also observe that the identification made in Eq. (^ seems to require a large number of particles 
A'", in order to make physical sense (see Ref . ) . 



3 Discussion 



The purpose of this section is to put forward possible problems/open questions in our approach, 
in order to stimulate further research. 

1. First of all, we have tried to look into the results of Refs.||l|, ^, ^ from a different perspective, 
and we have found that (for the purpose of this work) the assumptions of those papers could 
be replaced by others. Indeed, rather than saying that our systems are Anosov-like (or axiom 
A-like), we merely use one representation of the average currents, Eq.(|^, which would be 
correct for all smooth functions if the systems were really Anosov, and we postulate that the 
currents associated with UPOs are well behaved. This is similar to equilibrium statistical 
mechanics where rather than specifying the expected chaotic properties of the given systems, 
the form of the ensembles are postulated, and calculations are developed afterwards with such 
ensembles. One advantage of that might be that we do not really have to know which properties 
of the dynamics of particle systems make them "Anosov-like". In particular, smoothness of 
the dynamics, hyperbolicity etc. may not be strictly necessary, as several examples seem to 
indicate, especially if only selected functions of phase are considered. For instance, the Lorentz 
gas has singularities (although it is strongly hyperbolic); while the Ulam map is not hyperbolic. 
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Yet, POE have performed correctly when apphed to such systems (see Ref. p7[| for the Ulam 
map). See also Ref. [28|, in which intermittent diffusion is treated in terms of cycle expansions. 

2. Another advantage of our approach is precisely the fact that the special representation of the 
stationary state which we use —that based on orbital measures— can actually and succesfully 
be implemented in the simplest particle systems (various periodic Lorentz gas models), at 
least in the calculations of several quantities |15|, |16|, |l^, ^l|. For the same reason, our 
predictions formulated in terms of POE are amenable to direct tests via computer simulations 
of such systems. 

3. In all approaches, ours as wel as those of Refs.||l], ^, the necessity of having a large 
number of particles seems to be absent. Classical proofs of the Onsager relations always relied 
upon the assumption that the system is large (made of ~ 10^^ particles), and that it can be 
split into small "local" ones; basic variables ai, a^, with n <C A were associated with these 
small local subsystems (n ~ 10^'^ say). The origin of this apparent paradox is the assumed 
equality of the average phase space contraction and macroscopic entropy production rates, 
which only seems to hold for large N, Ref.|26|. 

4. In the paper we have avoided the problem of dealing with the singular points of the 
dynamics of systems of hard particles, saying that they would only get vanishing weights in 
the POE Eq.(0). In the Appendix, we give an argument to support this point of view. If that 
is correct, and orbits with grazing collisions effectively do not spoil the POE-representation of 
the stationary distributions, because they do not contribute to Eq.(|^, then we have evidenced 
one possible mechanism through which the techniques devised for axiom-A systems remain 
valid in more general settings, e.g. those of Refs.[@, 15, O, 17, 18, 21]. 



Appendix 

Here we argue that in case of hard core interactions, UPOs with grazing collisions can be 
neglected in the expansion Eq.(|^, so that they do not endanger the validity of Assumptions 1. 
and 2. First of all, the largest Lyapunov exponent of one such orbit would appear to be very 
large (effectively infinite) in numerical simulations; hence a very small (effectively vanishing) 
weight would be assigned to such orbits (see e.g. Ref.[^). Indeed, consider the trajectories 
of two particles which are going to experience a grazing collision. If a perturbation moves the 
trajectories of the two particles a little away from the colliding path, there will be no collision 
at all, while there will be a non grazing collision if the perturbation moves the trajectories 
closer to each other. 

Moreover, for the Lorentz gas at equilibrium, these ideas can be supported by an analytical 
argument. The Lyapunov exponents, in this case, are obtained from the logarithms of the 
eigenvalues of products of matrices of the following form:t3 

n 

if n is the number of collisions, where 

\ J cos tpj / 

^"This follows from a simple calculation. See, e.g., Ref. |30) . 



(29) 



(30) 
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Here, Ij is the distance travelled by the moving particle between the collisions j and j + 1, 
and ipj is the collision angle at the jth collision, where ipj = means head on collision, while 
ipj = 7r/2 corresponds to a grazing collision. Because ^ and cos ^pj > for all j, 

every A4j is defined and consists of negative entries if 'i/'j £ [0,vr/2) (hard collisions), while it 
is not defined for grazing collisions. 

Let us compute A4 for an orbit with n collisions, and ipj < 7r/2 for all j = 1, ...,n, and 
denote by 

M(-''^ = l[M,^{-l)-'(^l ^) , (31) 

the product of the first (n — 1) blocks. The terms a, b, c and d are then finite, positive numbers. 
Now, multiply A^^""^-* by Ain, the matrix of the last free flight and collision, in order to obtain 
M: 

cos i>n 

The eigenvalues of M. obey 



( -l\n \ ■ COSl/'n COSl/'n \ /ooN 



COS Ipn 



K + ^s = Tr{M) = {-lT [a + d + hln + -{c + dln)] , (34) 

V cos Ipn J 

and 

A^A^ = det(X) = ±1 . (35) 

Hence, the determinant remains bounded, while the magnitude of the trace grows without 
limits if ipn '?i"/2, which implies that one of the eigenvalues becomes very large, while the 
other becomes very small. The result is that the weight A~^ in Eq.(|^ gets smaller and smaller, 
and converges to zero for UPOs with ^„ closer and closer to 7r/2. 

This argument is not yet a full proof that orbits with grazing collisions can be neglected 
in the POE, since we have not explained in which sense one such orbit could be considered 
as the limit of a sequence of orbits without grazing collisions. However, the argument shows 
how a collision which is sufficiently close to grazing contributes to make the weight of the 
corresponding orbit small and, keeping fixed the other quantities, the weight is the smaller, 
the closer to grazing the collision is. We conclude noting that the problem posed by trajectories 
with grazing collisions is present in all approaches based on dynamical weights. In this respect, 
one of the nice features of the approach based on periodic orbits is that each UPO, as a whole, 
covers only a finite length, hence the relevant stable and unstable manifolds, and the UPO's 
contribution to the POE can be computed in finitely many steps. This is why the previous 
analysis for the Lorentz gas can be carried out. 
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